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Abstract 

Although Planck data supports the standard ACDM model, it still allows for the 
presence of Dark Radiation corresponding up to about half an extra standard neu- 
trino species. We propose a scenario for obtaining a fractional "effective neutrino 
species" from a thermally produced particle which decays into a much lighter stable 
relic plus standard fermions. At lifetimes much longer than ~ 1 sec, both the relic 
particles and the non-thermal neutrino component contribute to Dark Radiation. 
By increasing the stable-to-unstable particle mass ratio, the relic particle no longer 
acts as Dark Radiation but instead becomes a candidate for Warm Dark Matter 
with mass 0(1 keV - 100 GeV). In both cases it is possible to address the lithium 
problem. 
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1 Introduction 



The recently announced first cosmological Planck results [1] herald a new era in cosmology 
in which the standard ACDM model can be tested to high precision. For example, Planck 
results [1] combined with WMAP polarisation and other CMB data measure the excess 
of radiation at recombination to be, in units of neutrino species, 

AT e p ff lanck = 3.36 ± 0.34 , (1) 

assuming a standard value of the primordial helium abundance. This result supports the 
standard ACDM model while not excluding the possibility of an extra radiation compo- 
nent, so-called Dark Radiation (DR), beyond the standard one. Indeed, when the Planck 
data are combined with the Hubble constant H measurement from astrophysical data 
sets, in particular from the Hubble Space Telescope (HST), the best-fit value increases to 

^Planck = 3 q 2 ± Q 25 ^ ( 2 ) 

which amounts to a 2.3 a signal for DR. This shows that scenarios beyond the traditional 
(one fully thermalised) sterile neutrino hypothesis, in which a fractional "effective neutrino 
species" can emerge, are not excluded by Planck results and may even be mildly favoured 
by some data sets. 1 

In this paper, motivated by the above considerations, we shall propose a scenario 
for obtaining a fractional "effective neutrino species" from a thermally produced particle 
which decays into a much lighter stable relic plus a non-thermal active neutrino com- 
ponent. On the other hand, by increasing the stable-to-unstable particle mass ratio to 
0(0.1), the stable relic no longer acts as Dark Radiation but instead becomes a candidate 
for Warm Dark Matter (WDM). Thus our scenario is flexible enough to account for either 
DR or WDM (but not both at the same time). Interestingly, in both cases it is possible 
to address the lithium problem. However, before discussing details of our scenario, it 
is worth recalling some general constraints on new physics beyond the standard ACDM 
model. Although well known, they are worth recalling at this point since they provide 
important constraints on our scenario. 

Big Bang Nucleosynthesis (BBN) is the most traditional cosmological probe of new 
physics [6] . Non-standard BBN effects have been extensively studied within scenarios pro- 
ducing modifications of the neutrino content compared to the Standard Model (SM) [7] 

1 For a recent discussion (in the light of Planck) on models able to yield fractional "effective neutrino 
species", see [2]. Notice that also in active-sterile neutrino oscillations, beyond a full sterile neutrino 
thermalisation, fractional "effective neutrino species" can be obtained [3] and in this case the Planck 
data impose stringent constraints on the mixing parameters that seem to indicate a strong tension with 
the short-baseline hints [4]. A traditional solution, though difficult to justify is to assume large initial 
lcpton asymmetries suppressing the mixing [5]. 
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such as a massive (m v = 0(10 Me V)) decaying ordinary neutrino, now excluded by neu- 
trino oscillation experiments [8], or active-sterile neutrino oscillations [3]. On the other 
hand the inclusion of ordinary neutrino oscillations does not produce any modification 
to the standard scenario, in particular the SM value iVJg — 3.046 does not change [9]. 
Scenarios where massive metastable (r > 100 sec) particle decay products modify the 
abundances after nucleosynthesis have also been extensively investigated [10]. 

However, the possibility to test non-standard BBN effects is greatly limited by system- 
atic uncertainties in the determination of the primordial nuclear abundances [11]. The 
discovery of the acoustic peaks in the power spectrum of CMB anisotropies has opened 
new opportunities to probe non-standard effects with much lower systematic uncertainties. 
First of all, from CMB, it has been possible to measure with great accuracy and precision 
the baryon-to-photon number ratio t]b that, if assumed to be constant between BBN and 
recombination time, allows to make firm predictions on the Standard BBN (SBBN) values 
of the primordial light element abundances to be compared with the measured values. 

Moreover, from the observed acoustic peaks it is possible to constrain the presence 
of DR at the recombination time, the hot Dark Matter (DM) contribution, and even 
the primordial helium-4 abundance Y p . Recent data from CMB observations have hinted 
to non-vanishing DR [12] and this has triggered quite an intense investigation on the 
possible sources, ranging from sterile neutrinos [13, 14] over modifications of the neutrino 
temperature [15] to exotic relativistic species [16]. 

Non-standard effects from the decays of long-lived massive particles have been quite 
extensively investigated in general [17], in the case of inert product particles [18], and 
in the case of electromagnetically interacting product particles, when they can alter the 
primordial abundances after nucleosynthesis [10]. In this case these have been advocated 
to reconcile a tension between the observed lithium abundance that is about three times 
lower than the value predicted by SBBN, the so-called lithium problem [19]. 

Returning to our scenario, this needs to be considered carefully since it could poten- 
tially produce non-standard effects in a non-trivial way. For instance, the decay of a 
heavy thermally produced particle species into a new weakly coupled lighter stable relic 
plus non-thermal neutrinos, after the freeze-out of the neutron-to-proton abundance ratio 
and the neutrino decoupling, could alter the lithium abundance without affecting signif- 
icantly the helium abundance. Concerning the details of our scenario, it is based on the 
particle physics model introduced in Ref. [20], where two Majorana fermions xi,2 with 
masses M 2 > Mi are added to the SM and both fermions can couple to the Z-boson, 
but only with vertices suppressed by factors e 1}2 - However the cosmological implications 
of this model that we consider here are completely new. In particular the possibility for 
DR has not previously been considered, and the mechanism for WDM production here is 
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quite different from the previous case where the xi was produced in thermal equilibrium. 
In the present case, only the X2 is produced in thermal equilibrium and decays with a 
lifetime r as X2 — > Xi + / + /> where / is any SM fermion, which (depending on the 
parameters) allows for either DR in Xi an d v or WDM in Xi- DR emerges if the light 
stable relic xi is much lighter than X2, while WDM is obtained by increasing the mass 
ratio to 0(0.1). 

The plan of the paper is the following. In Section 2 we discuss the basic features of 
the model. In Section 3 we calculate the DR contribution and show that this can explain 
the required value of AN cS = N cS — 3.046 from the combination of Planck and HST data. 
In Section 4 we discuss the alternative scenario where the model provides an explanation 
for DM while the amount of DR would be negligible. In Section 5 we discuss how, in 
both cases, the lithium problem can be addressed, though for different ranges of values of 
the mass of the decaying particle. Finally, in Section 6, we draw the conclusions. In the 
Appendix we provide some technical details on the decay rate. 

2 The model 

Our setting extends the considerations of Ref. [20]. We assume, in addition to the SM, two 
Majorana fermions xi,2 with masses Mi )2 . These fields are mainly SM-singlets, but have 
interactions with the Z-boson which are suppressed by factors ei )2 and 5. The interaction 
Lagrangians are given by 

• Z-Xi-Xi (* = !> 2 ) : 

Ca = gZ^x~a^5Xi, (3) 

• Z-xi-X2- 

£12 = ge^xlYiA^L + B X P R ) X 2Z^ + /i.e., (4) 

(Pr,l = (1 ±75)/2) with g ~ 0.653 being the generic SU(2) gauge coupling and (A X ,B X ) 
with A 2 X + B X = 2 being constants which parametrize the structure of the coupling. While 
at this stage these Lagrangians are only a postulate, similar settings are known to exist 
for example in the E^SSM [21] or in certain Left-Right symmetric models [22]. 

The most important point is that such a setting admits a decay X2 — > Xi + Z*, Z* — > 
//, where / is any SM-fermion. As long as the masses M 12 are suitably chosen, this in 
effect amounts to a transition X2 — > Xi + f + f- 2 Then, we can make the simple observation 

2 As illustrated in Rcf. [20], the other possible mode \i ~^ XiW + W~ does not make much of a 
difference, since it is either kincmatically forbidden or just a small perturbation. The mode \i ~^ ^Xi is 
even further suppressed. For simplicity, we neglect both of them. 
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that this reaction in particular allows for a decay involving light neutrinos 

X2 ->• Xi + v + V, if M 2 > Mi. (5) 

If we furthermore make the assumption that e 2 is large enough to keep xi in thermal 
equilibrium in the early Universe, while e\ is sufficiently small that this does not happen 
with xi, we can have the following phenomenologically interesting cases: 

1. Ordinary neutrinos and Xi's as Dark Radiation: 

In a certain parameter range we produce DR. Remarkably a non-negligible fraction 
of this DR consists of ordinary neutrinos, which might have further interesting 
implications. In other words, on top of the thermal neutrino component in the 
early Universe, the decays produce also a non-thermal contribution, similarly to a 
scenario discussed in [23]. 

2. Xi as Warm Dark Matter: 

When a sufficiently large abundance of X2-particles freezes out, they will all decay 
and each decay produces exactly one Xi- For the right combination of masses and 
couplings, xi could play the role of WDM if it is not too hot. 

We will now analyse both situations from a phenomenological point of view. It will 
turn out that indeed both points can be fulfilled in certain regions of the parameter space, 
however, they do not work out simultaneously. In other words, in the parameter regions 
where we obtain the correct WDM abundance we have practically no DR, and in the 
regions where we get reasonable amounts of DR the Xi' s also contribute to DR but not 
to WDM. 

3 The Dark Radiation scenario 

In this section we calculate AiV e gf, defined as the value at recombination, within our 
model, to be compared with the value in Eq. (2), found by combining Planck data with 
the Hubble constant measurements. 

The total energy density pr in the radiation component receives a contribution from 
standard particles and from the Xi's- In our case, if we restrict ourselves to the case where 
the helium abundance is standard, corresponding to setting r to values much longer than 
the neutron-to-proton number ratio freeze-out time tf T ~ 1 sec, the neutrino contribution 
can be unambiguously split into a standard thermal component p* h and into a non-thermal 
component p° th resulting from the Xi decays. The X2's are non-relativistic at the time 
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of their decays and they have completely disappeared at the recombination time t rec . 
Therefore, the radiation contribution can be written as 

PR (T)=g R (T)^T\ (6) 

where the number of radiative degrees of freedom can be expressed as the sum of a SM 
component and of a non-standard component given by the non-thermal neutrinos and by 
the xi contribution so that g R {T) = g% M (T) + gf h (T) + g xl (T). The SM contribution is 



9 S R M (T) = 2 + 7 - 



9e4T) + 2[^) iV c s ff M (T) 



(7) 



At temperatures T <^ m e the contribution from e ± vanishes, the neutrino-to-photon 
temperature ratio saturates to its asymptotical value T u /T = (4/11) 1 / 3 ~ 0.715, and the 
effective number of neutrino species freezes to the value TVjjj 1 ~ 3.046, differing from 3 
since the thermal neutrino component is actually very slightly heated by e ± annihilations. 

The non standard component can be analogously parametrised in terms of the extra 
number of effective neutrino species 
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rnthm , mi 4/ny i2o (ny Pdr (t) 

AN cS {T < m e ) = [g v (T) + g Xl (T)\ - I — J =^ I jl y4 , (8) 

where we defined /?dr = Pxi + P° th - The energy densities of the non-thermal neutrino and 
Xi components obey very simple fluid equations, 

d -^p- - j(p x ,R 3 ) - (pTR 3 )h (9) 

and 

^7^ = ^(PX^ 3 )-(PX^ 3 )^ (10) 

where b u = 2BR^/3 and b Xl = 1/3 are, respectively, the averaged fractions of energy into 
neutrinos and Xi an d BRj, is the branching ratio of Xi decays into neutrinos. 

Notice that we are assuming M 2 3> Mi so that the xi can be treated as ultrarelativis- 
tic at the production. However, in order for the DR contribution not to be negligible, 
the Xi have to be necessarily ultrarelativistic not only at the production but even until 
recombination since otherwise they would over-contribute to the DM energy density. 3 For 
this reason the xi's contribute to AA^fr until recombination. 



3 This can be seen in a qualitative way imposing that the average momentum of Xi & a t icq (the 
matter-radiation equality time) p cq » M 1; leading to the condition Mi/M 2 <C 10~ 6 ^/r/sec. 
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Assuming radiation dominance until equality, 4 a solution of the differential equations 
is quite straightforwardly found, 



PBnR 3 = b BR m X2 N l X2 J dt'^— a ^ = b BR m X2 N l X2 ^^m, (H) 
having defined &dr = K + b Xl and introduced 




where the error function is defined as erf(x) = (2/y / 7r) J^e 2 dz, such that simply 
f (t) ^ 1. From Eq. (8) we can then calculate 



-0.47 fed(t )6 DR1 ^^-«t), (13) 

where we introduced the dilution factor d(t) = N^/Nj(t) and the %2 relic abundance iV"£ 2 
that got frozen at some freeze-out time tf. Notice that in the numerical expression we 
have approximated g]( 4 ~ (fi^ 1 ) 1 ^ 4 , neglecting a small correction (~ 2%) from the DR 
itself. We have also used a normalisation of the scale factor such that the abundance of 
X2 in ultrarelativistic thermal equilibrium is just 1 (i.e. N^{T ^> M 2 ) = 1). 

The asymptotic value of Eq. (13) is then simply obtained by taking £ = 1 so that 5 



M 2 N X2 , T 



AN eS (t > r) ~ AiV cfr ~ 0.47 ^ d &DR -r^rf */ — . (15) 

MeV V sec 



4 Notc that, in the region where we obtain a suitable value of AN c g, \i can however never dominate 
the energy density of the Universe, as otherwise it would produce by far too much DR. 

5 This equation agrees with the result obtained in Ref. [24] for decays just into inert particles, which 
would be formally recovered setting &dr = 1 and taking into account that in our normalisation 

rf _ Y L 9S0 8^ 4 
d g X2 135 C(3) 

where Y£ 2 = n X2 /s, s is the entropy density, gso — 3.91 is the entropy number of degrees of freedom 
at the present time, and do is the dilution factor at the present time. Notice that our model does also 
include this (more traditional) scenario that would be recovered for lifetimes r -c t^ cc < tf T , where t^ ec 
is the time when standard neutrinos decouple. In this case one would then have simply p" th = 0. For 
lifetimes t ~ t^ cc , tf T neutrinos would be in the decoupling stage and, therefore, this range would require a 
more complicated kinetic analysis both for the calculation of the contribution to AN c g from non-thermal 
neutrinos and, as we will point out in Section 5, also for the primordial helium abundance. 
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Figure 1: Dark Radiation from \ 2 — > Xi + / + / a ^ recombination time t rec for a lifetime 
r = 10 sec. 



In Fig. 1 we show both the total contribution AN c g to DR (left panel) and just the 
neutrino contribution (6y/6 DR ) AA^ efT (right panel), for different values of the mass M 2 
as functions of In the plots, we have ei < 10~ 4 , such that Xi never enters thermal 
equilibrium and is only produced by X2-decays. We have also chosen r = 10 sec for the 
plots, which can always be achieved by varying the remaining parameter 5, cf. Appendix A 
for more details. 

The blue/light gray regions mark the la regions from the Planck plus other CMB data 
(dashed horizontal lines), as well as the corresponding regions when the if - measuremen t s 
are also taken into account (solid horizontal lines), cf. Ref. [1]. The dark gray region on 
the right is excluded by the invisible decay width of the zT-boson, which constrains e 2 to 
be smaller than about 0.23 for e\ <C e 2 , in case that M 2 < Mz/2 [20]. As can be seen, 
for large enough masses M 2 we can always find a region of e 2 where we have a significant 
amount of DR while still satisfying the bounds. This contribution becomes large for small 
enough e 2 , since this region corresponds to a very early freeze-out of xi an d hence to an 
unsuppressed abundance. 

Fig. 1 might suggest a laboratory detection of the additional amount of ordinary 
neutrinos via neutrino capture on tritium or rhenium [25] or via (potentially resonantly 
assisted) modified electron capture on holmium [26, 27]. Unfortunately, the rates turn out 
to be very small, due to the additional neutrino component in spite of the higher energy 
still being much less than the ordinary cosmic neutrino background. Even resonance 
enhancements [26] are not powerful enough to alter this conclusion. 
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4 Stable relic as Warm Dark Matter 



In Ref. [20] it was shown that a thermal production of stable Xi's would typically overdose 
the Universe for small masses Ml, but this can be cured by entropy production from \2 
decays to yield the keVin/GeVin scenario for DM. We take a different route here by 
assuming that e\ is small enough (< 1CT 4 ) to prevent Xi from ever being in thermal 
equilibrium. We assume that the xi s are entirely produced from xi decays, so that we 
have a new way to produce a suitable DM candidate. Similar settings can be found in the 
literature, see e.g. Refs. [28]. Since our scenario is more constrained than some of them, 
we have an important bound originating from the invisible Z-boson decay. 

As had been shown in Ref. [20], one can thermally produce a non-relativistic abundance 
of %2- 6 The \2 abundance is, as for a generic non-relativistic species, proportional to the 
mass M 2 . However, xi is unstable and decays with a lifetime r ^> tf T . Since each xi 
decay produces exactly one xu this allows to simply translate the xi abundance into the 
abundance of the DM candidate xi ( & t the present time): 

n DM h 2 = n xl h 2 = ^n X2 h 2 , (16) 

where Q X2 h 2 is the final abundance of X2 if it was the DM. By Eq. (16), we can always 
correct an overabundance in X2 such that the final DM abundance in Xi hits the observed 
value, tt DM h 2 = 0.1196 ± 0.0031 [1]. 

Note that xi is produced in the early Universe at times t ~ r, but the expansion 
redshifts its momentum and, for this reason, even if it is ultra-relativistic at the produc- 
tion, it can still be slowed down by the cosmic expansion and be a viable DM candidate. 
We take this into account by calculating the free-streaming scale of Xi- This calculation 
can be found in textbooks (see, e.g., Ref. [29]), and it amounts to calculating the mean 
distance which the particle would travel if it was not trapped gravitationally. Technically, 
one has to evaluate: 

A FS (t)= fdt V -f (17) 

where v(t) is the velocity of xi and a(t) is the scale factor. Using elementary kinematics 
and the approximation of radiation-domination until equality, we obtain 

W() K o, Mpc ^ {^l) h. (VA + VTTA) , (18) 

6 In principle, \i could also be relativistic at the production, but this would impose further complica- 
tions and most probably not improve the result, since very hot particles would be produced. 
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where 



( 



M 2 /M] 
10 2 




2 



-1 



A ~ 0.181 x 10 ; 



,8 



10 sec 



T 



) 



(19) 



In order for the smallest structures in the Universe not to be erased, we need Aps ^ 
0.1 Mpc [30], and hence we must be in the region 



(of course Mi/M 2 < 1). A free-streaming scale of ~ 0.1 Mpc corresponds to WDM, 
while values much below would correspond to cold DM. Therefore, for minimum values 
t ~ 10 sec one has already a marginal allowed region for cold DM (M 1 /M 2 > 0.1), that 
tends to disappear for longer lifetimes. For this reason in this setup DM is typically warm 
rather than cold. 7 

We have calculated the abundance of \ 2 along the lines of Ref . [20] , but with a more 
precise version of the decay width, cf. Appendix. The results can be seen in Fig. 2, where 
we have plotted the lines of correct abundance for different values of M 2 as functions of the 
suppression parameter e 2 . Indeed, we are hit by the bound from Z-decay such that both 
M 2 and M\ must be relatively large, which is very different from the scenario presented in 
Ref. [20]. A particularly interesting point is that even a relatively heavy xi of 100 GeV, 
or so, could be a warm species produced at temperatures around 1 MeV. The intuitive 
reason is that the energy that it obtains was in some sense "stored" at a relatively low 
temperature inside the non-relativistic %2's before they started decaying. Of course notice 
that the %2's do nevertheless decay much before matter-radiation equality and, therefore, 
they do not act as an additional DM component. 

At this stage one could wonder whether the DM solution is compatible with having 
a sizeable DR contribution (calculated in the previous section), certainly an intriguing 
possibility. In this case the xi s would be already non-relativistic before matter-radiation 
equality and, therefore, only non-thermal neutrinos would contribute to DR. Unfortu- 
nately, as already anticipated at the end of Section 2, the model cannot lead to a sizeable 
amount of DR in the same region of parameter space where one has the correct Aps, the 
region marked in pink/light gray in Fig. 2. This is easy to understand: a large value 
of e 2 keeps the %2's in equilibrium for a long time, such that their freeze-out abundance 
iV£ 2 is suppressed. While they can nevertheless produce a significant non-relativistic en- 
ergy density in xi's, due to the large mass Mi, the mere amount of neutrinos produced 
by the reaction quoted in Eq. (5) is not enough to contribute significantly to DR given 

7 Noticc, however, that in the case r < 10 sec, that we are not considering (cf. footnote 2), the cold 
DM region would become less marginal. 
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Figure 2: The masses Mi required to obtain the correct WDM abundance after the 
decay of all (with r = 10 sec), displayed for different masses M2 as functions of the 
suppression parameter €2- The invisible Z-decay width requires e 2 < 0.23 for M2 < Mz/2. 
Note that we get warm DM in the marked region, while CDM is only marginally possible 
due to e 2 < 1. 

the observed value of Q-Duh 2 . This means that a sizeable value of AiV e g ~ O(0.1) is 
incompatible with the DM solution. 

This can be easily understood also on quantitative grounds. Indeed if one plugs the 
DM condition Eq. (16) into the Eq. (15), one finds the following relation linking AiV e ff to 



One can then immediately see that, in order for the condition Eq. (20) to be satisfied, 
one needs AN e n < 0.001, saying that the DM scenario predicts, within conceivable ex- 
perimental precision, a vanishing DR contribution. 

Notice that there are two further restrictions on r: (i) it cannot be smaller than the 
time of the freeze-out of \i an d also (ii) the condition 8 < 1 enforces r not to be smaller 
than a certain minimum value depending (mainly) on M 2 (see the Appendix). Finally, 
on top of the further reaching cosmological aspects of our setting, we would like to add 
that one could even think of discovering our model at LHC, by directly producing \2~ 
pairs at high enough energies. We will postpone a more detailed investigation of all these 
interesting aspects to future work. 



Mi/M^ 



2- 




(21) 
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5 Addressing the lithium problem 



The observation of acoustic peaks in the power spectrum of the CMB anisotropies gives a 
very precise measurement of the baryon-to-photon number ratio f]^ MB = (6.030 ±0.075) x 
10~ 10 [1]. This can be used to derive quite precise predictions of the primordial nuclear 
abundances of light elements within SBBN that can be compared with the observed ones 
in astronomical environments. 

The primordial nuclear abundances that can be used as cosmological probes are 
the helium-4 abundance Y p , the deuterium abundance D/H and the lithium abundance 
Li/H [31]. In the first case, with an (careful and significant) enlargement of the systematic 
uncertainties one finds, from clouds of ionized hydrogen (Hll regions) in dwarf galaxies, 

Y p = 0.249 ± 0.009 (95% C.L.) , (22) 
in agreement with the SBBN prediction [32] 

^sbbn^cmB) _ a2466 + Q1 1^(^/5) ~ 0.2474 , (23) 

where 7710 = lO 10 ^- 

The primordial deuterium abundance, measured from high-redshift, low metallicity 
quasar absorption systems, is found to be 

(D/H) = (2.82 ± 0.21) x 10~ 5 (95% C.L.) , (24) 

also in agreement with the inferred SBBN value 

(D/H) SBBN (^ MB ) ~ 3. 6 x io-5 foo/5)- 1 - 6 ~ 2.7 x lO" 5 . (25) 

On the other hand a stellar determination of the lithium abundance gives (Li/H) = 
(1.7 ± 0.06 ± 0.44) x 10~ 10 , about twice lower than the SBBN prediction (Li/H)p BBN ~ 
(3 — 5) x 10~ 10 , the mentioned lithium problem that, barring still unknown conventional 
astrophysics mechanisms, provides a potential evidence of non-standard effects. 

It has been shown [19] that the decays (into hadrons and/or muons) of long-lived 
particles with a lifetime r = 0(1O 4 sec) and abundances 0(1) — 0(100) times the baryon 
abundance, can reduce the lithium abundance by 0(1) factors, thus solving the lithium 
problem, without spoiling at the same time the agreement of the deuterium and helium 
abundances. 

These features can be easily fulfilled by our \i decays. If we impose the necessary 
conditions on the lifetime r and on the relic abundance N* , we can show that we obtain 
physical solutions both in the DR and in the DM scenario. 
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Let us start with the DR scenario. First of all, the relic abundance iV£ 2 is related to the 
X2-to-baryon number ratio i] X2 = n X2 /n B at r ~ 10 4 sec by the relation N %2 ~ 10~ 8 i] X2 . 
Therefore, plugging this relation into Eq. (15) with r ~ 10 4 sec and imposing the Planck 
value AN cS ~ 0.5 (cf. Eq. (2)), one finds easily M 2 = O(10 GeV) - 0(1 TeV). 

Let us now consider the WDM scenario. Repeating the same calculation with the 
condition AN„ < 0.001 found in the WDM scenario, one obtains M 2 = 0(100 MeV) - 
0(1 GeV), having taken into account that necessarily \2 should be heavier than muons 
for the solution of the lithium problem to be viable. 

These results show that our model can nicely connect two independent cosmological 
puzzles, respectively, either DR or WDM and the lithium problem. 

In the case of the DR scenario there is also another interesting consideration to be 
done. For interesting values of the lifetime r = 0(1O 4 sec), the model still predicts a SBBN 
value of Y p since AiV e ff(tf r ) would be still vanishing. Therefore, if future data supported 
values AiVefr — 0.5, then it would be interesting to be able to test AN e g(tb) "C 0.5, 
since in this way the model could be distinguished from a more traditional case where 
DR is generated prior to t fr and then constant. The SBBN primordial helium abundance 
prediction would be modified by a quantity [32] AY P ~ 0.0137 AiV cff . From Eq. (22) one 
then finds the bound AiV eff (tf r ) < 0.75, that is still currently not stringent enough to 
provide an indication for such a scenario. It is interesting that CMB observations are also 
able to provide a measurement of Y p , though currently the error is even larger than those 
from astronomical data sets Eq. (22). If future measurements will be able to constraint 
AY p <C 0.006 this could provide then another interesting piece of information. 

Analogous considerations can be done for deuterium. In this case this is sensitive 
to AiVeff(t n uc), where t nuc ~ 365 sec is the time of nucleosynthesis. For life times r ~ 
0(1O 4 sec) one would expect AiV eff (t nuc ) <C AiVefr, though in this case one could have some 
small contribution if t ~ 0(1O 3 sec). Then a non- vanishing AN eS (t nuc ) would modify the 
SBBN into [32] 

(D/H) ~ (D/H) SBBN (^) [1 + 0.135 AN cS (t QUC )] - 8 . (26) 

Therefore, comparing the observed value Eq. (24) with the SBBN prediction Eq. (25), 
one finds a constraint AiVeff(t nuc ) < 1-0 that is also too large to draw a conclusion on 
a possible difference between AiV eff (t nuc ) and AiV cff . However, a future improvement in 
the determination of D /H might make that possible if at the same time a non- vanishing 
AiV c ff ~ 0.5 should be established. 

These considerations show how, in the case that future data should give an evidence 
for a non-vanishing AiV e ff ~ 0.5, then a combined analysis with light element primor- 
dial abundances can be potentially be used to highlight a dynamical evolution of DR as 
predicted by our model. From this point of view the current lithium problem would be 
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interpreted as an effect of the X2-decays. 

6 Conclusions 

We have seen that while Planck results support the standard ACDM model, they still 
allow for an amount of DR corresponding to about half an extra neutrino species, even 
favoured if data from astrophysical data sets (e.g. the HST) are taken into account. In 
particular, one could argue in support of cosmological non-standard effects that seem to 
require a dynamical mechanism for the production of Dark Radiation at times after the 
freeze-out of the neutron-to-proton number ratio. 

Motivated by this situation, we have proposed a scenario for obtaining a fractional 
"effective neutrino species" from a thermally produced particle which decays, much after 
neutrino decoupling, into a much lighter stable relic plus a non-thermal active neutrino 
component. 

Our scenario is based on two Majorana fermions, the heavier of which can decay into 
the lighter one plus a pair of ordinary fermions. If the lifetime is much longer than 
the neutrino decoupling time, the DR consists of two components, the lighter Majorana 
fermion and a non-thermal ordinary neutrino component and there would be an interest- 
ing, potentially testable, dynamical evolution of DR, with different values of AN eS (t) at 
the different observationally relevant times. For this reason we have limited our analysis 
to this range of lifetimes (r»l sec). 

We found that a sizeable amount of DR, within the reach of future experimental 
investigations, is produced when the stable relic is orders of magnitude lighter than the 
decaying particles. On the other hand , by increasing the mass ratio M1/M2 to 0(0.1), the 
stable relic no longer acts as Dark Radiation but instead becomes a candidate for WDM. 
In this way we obtained two mutually exclusive scenarios, one predicting a sizeable amount 
DR and one explaining DM. 

Interestingly both scenarios, for r ~ (9(10 4 sec), are compatible with a solution of the 
lithium problem due to a partial disintegration of the synthesised value operated by the 
decay products if the mass of the decaying thermal relic particle that should be within 
the interesting range M 2 = O(10 GeV) — 0(1 TeV). Hence in both cases it is potentially 
possible to resolve the lithium problem. In the case of DR, this would also predict a 
SBBN and close-to-standard value respectively for the primordial helium and deuterium 
abundances (i.e. vanishing AN c g(t) at tf r and AN e g(t) smaller than AN e g at t nuc ). 

In conclusion, we have seen that a relatively minimal extension of the Standard Model 
of particle physics is capable of predicting a few new signals in cosmology beyond the 
standard ACDM model. The proposed scenario seems to provide quite a flexible frame- 
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work which could account for either DR or WDM. In the light of the first cosmological 
Planck data we have seen that there is a hint for DR and it is possible that a positive 
signal could be observed in future results. We eagerly look forward to the next results 
from Planck in the new era of precision cosmology in which scenarios such as the one 
presented here will be fully tested. 
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Appendix: Decay width for \i ~ * Xi + v + v 

The decay width looks easiest for the case of the final state fermion / being a neutrino, 
since then it is effectively massless. From Eq. (4), one can derive 8 

r(x 2 xi**) = ^f^f 1 • A l ■ t 2 ^ - H x l 3 ], (A-l) 

M 2 / -. r- _ M 2 

where 1 5 = J dt g(M 2 , M u 0, t) ■ t\ 1 + and X 3 = / dt g(M 2 , M u 0, t) ■ \t 2 . 



t=0 ,21 t=Q 



Furthermore, we have defined A v = 9y 2 ^ A , F x = A 2 X + B 2 ,, and H x = 2A X B X , for an 
interaction Lagrangian C = g^xil^{.A x P L + B x P R )x2Z^ + h.c. Finally, the integrand 
functions g(M 2 , Mi, to/, t) and $(M 2 ,M 1 ,t) are defined as 



>t 2 - 4m 2 f V$ 

9 - (t2 V _ M|)2 + r|MF * - K - (* - m M ~ (t + MJ] . (A-2) 
These are the expressions we have used to obtain our numerical results. 
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